The shape of multidimensional gravity 
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In the case of one extra dimension, well known Newton's potential tp{rz) = —Gjim/rz is 
generalized to compact and elegant formula <p{rz,£) = — (GNm/rz) sinh(27rr3/a)[cosh(27rr3/a) — 
cos(27r^/a)] _1 if four-dimensional space has topology R 3 x T. Here, r$ is magnitude of three- 
dimensional radius vector, £ is extra dimension and a is a period of a torus T. This formula is 
valid for full range of variables r% £ [0, +oo) and £ £ [0, a] and has known asymptotic behavior: 
ip ~ I/V3 for r3 >> a and tp ~ l/rf for r4 = \J r \ + £ 2 << a - Obtained formula is applied to an 
infinitesimally thin shell, a shell, a sphere and two spheres to show deviations from the newtonian 
expressions. Usually, these corrections are very small to observe at experiments. Nevertheless, in 
the case of spatial topology R 3 x T d , experimental data can provide us with a limitation on maximal 
number of extra dimensions. 
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Introduction From school years we know that newto- 
nian gravitational potential of a body with mass m has 
the form of ip(r) = —G^mlr where Gm is Newton's grav- 
itational constant and r = |r| is magnitude of a radius 
vector in three-dimensional space. This expression can be 
derived from Gauss's flux theorem in three-dimensional 
space. However, it is possible that spacetime might have 
a dimensionality of more than four and posses a rather 
complex topology. String theory and its supersymmet- 
ric generalizations - superstring and supergravity widely 
use this concept and give it a new foundation. There- 
fore, it is of interest to generalize the Newton's formula 
to the case of extra dimensions. Clearly, the result de- 
pends on topology of investigated models. We consider 
models where (D = 3 + <i)-dimensional spatial part of 
factorizable geometry is defined on a product manifold 
Mr) = R 3 x T d . R 3 describes three-dimensional fiat ex- 
ternal (our) space and T d is a torus which corresponds 
to ci-dimensional internal space. Let Vd be a volume of 
the internal space and b ~ V^ d is a characteristic size 
of extra dimensions. Then, Gauss's flux theorem leads 
to the following asymptotes for gravitational potential 
(see e.g. [l(): ip ~ l/r 3 for r 3 >> b and tp ~ l/r^ 
for r^ + d « b where r 3 and r^+d are magnitudes of ra- 
dius vectors in three-dimensional and (3 + e?)-dimensional 
spaces, respectively. In our letter, we generalize the New- 
ton's expression to multidimensional case with the shown 
above asymptotes and demonstrate that at least in the 
case d = 1 it has the form of compact and elegant formula 
via elementary functions. 

Multidimensional gravitational potentials In D- 
dimensional space, the Poisson equation reads 

&dVd = S d GvPd{td) , (1) 
where Sd — 2n D / 2 /T(D/2) is a total solid angle (square 



of (D — l)-dimensional sphere of a unit radius), Gx> is a 
gravitational constant in (V = D + l)-dimensional space- 
time and pd{yd) — m5(xi)S(x2)---S(xr}). In the case of 
topology M. D , Eq. ([T|) has the following solution: 
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D > 3. 



(2) 



This is the unique solution of Eq. (TTJ) which satisfies 
the boundary condition: lim <^c(rrj) = 0. Gravita- 

tional constant Gr> in (JlJ is normalized in such a way 
that the strength of gravitational field (acceleration of a 
test body) takes the form: —dtpD /dro = — GT>m/r D ~~ . 

If topology of space is R 3 x T d , then it is natural to 
impose periodic boundary conditions in the directions 
of the extra dimensions: </?r>(r 3 , £1, £2, ■ ■ ■ , £i, ■ ■ ■ , £d) = 
Vc( r 3,£i,6, ••.,& + a-i, ■■ - ,£d), i = l,...,d, where a» 
denotes a period in the direction of the extra dimension 
Then, Poisson equation has solution (cf. also with 
El): 



¥>i)(r3,£i, -,£d) = - 

x J2 ■•■ J2 exp 

ki— — 00 k^ — — oo 
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x cos l £1 ...cos £ d 

a\ I \ a d 



(3) 



To get this result we, first, use the formula i5(£i) = 



~t~ J2k=~oo cos (li~£*J an d> secon d, put the following re- 
lation between gravitational constants in four- and V- 
dimensional spacetimes: 



Sd 
S 3 



Gt> 



= Gn 



(4) 
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The letter relation provides correct limit when all a,; — > 0. 
In this limit zero modes fc, = give the main contribution 
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and we obtain ^n(r 3 , £1, £<j) ~* ~G N m/r 3 . Eq. (g]) 
was widely used in the concept of large extra dimensions 
which gives possibility to solve the hierarchy problem [l|, 
3]. It is also convenient to rewrite ^ via fundamental 
Planck scales: 



^ • M 2 

g M Pl(4) 



— M 2 + d 

— IV1 Pl(V) 



(5) 



where M Pi 



G 



-l/(2+d) 



(4) 



G 



-1/2 



1.2 • 10 ia GeV and M Pi m 



(T) 



are fundamental Planck scales in four and T> 
spacetime dimensions, respectively. 

In opposite limit when all a; — > +oo the sums in Eq. 
([3]) can be replaced by integrals. Using the standard 
integrals (e.g. from [J) and relation we can easily 
show that, for example, in particular cases d = 1, 2 we get 
desire result: <£ D (r 3 ,fi, ...,£<*)-> -G c m/[(£>-2) r^]. 

One extra dimension In the case of one extra dimen- 
sion d = 1 we can perform summation of series in Eq. 
©. To do it, we can apply the Abel-Plana formula or 
simply use the tables of series Q. As a result, we arrive 
at compact and nice expression: 



<^4(r 3 ,£) 
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(6) 



where r 3 € [0, +oo) and £ <E [0, a]. It is not difficult 
to verify that this formula has correct asymptotes when 
r 3 >> a and r± << a. Fig. 1 demonstrates the shape 
of this potential. Dimensionless variables rji = r 3 /a G 
[0, +oo) and r) 2 = £/a G [0, 1]. With respect to variable 
i]2, this potential has two minima at 772 = 0, 1 and one 
maximum at 772 = 1/2. We continue the graph to neg- 
ative values of 772 £ [—1, 1] to show in more detail the 
form of minimum at 772 = 0. The potential ([6]) is finite 
for any value of r 3 if £ ^ 0, a and goes to —00 as — l/r| 
if simultaneously r 3 — > and £ — > 0, a (see Fig. 2). 




FIG. 1: Graph of function 0(771,772) = v?4(r3,£)/(Gjvm/a) = 
— sinh(27r77i)/[77i(cosh(27T77i) — cos(27T772))]. 

Having at hand formula ([6]), we can apply it for cal- 
culation of some elementary physical problems and com- 
pare obtained results with famous newtonian expressions. 
For our calculations we shall use the case of £ = 0. It 
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FIG. 2: Section £ = of potential (©. Solid line is ^(771, 0) = 
— sinh(27T77i)/[?7i (cosh(27T77i) — 1)] which goes to — 1/771 (dot- 



ted line) for 771 
771 -v 0. 



+00 and to — 1/(-7T»7j.) (dashed line) for 



means that test bodies have the same coordinates in ex- 
tra dimensions. It takes place e.g. when test bodies 
are on the same brane. Also, to get numerical results 
we should define the size a of the extra dimension. If 
the standard model fields are not localized on the brane, 
then experiments give upper bound a < 10 _17 cm 
For this value of a the 5-dimensional fundamental Planck 
scale is M Pl ^ } > 10 n TeV (see Eq. ©). This value a 
can be greatly increased if we suppose that the standard 
model fields are localized on the brane. The gravita- 
tional inverse-square law experiments show that there is 
no deviations from three-dimensional Newton's law up 
to 2.18 • 10~ 2 cm [6]. Thus, we can take for this model 
a « 1(T 2 l(T 3 cm which results in M Pl{5) « 10 6 TeV. 
However, it is necessary to keep in mind that this case 
can be constrained by observations from supernova cool- 
ing (see e.g. 0). 

It is worth of noting that all formulas in this letter are 
applied to the Coulomb's law if electromagnetic field is 
not localized on the brane. It was not found deviations 
from three-dimensional Coulomb's law up to 10 _16 cm [5|. 
Therefore, for models with non-localized electromagnetic 
field we should take a < 10 _17 cm. 

Infinitesimally thin shell Let us consider first an in- 
finitesimally thin shell of mass m = 4irR 2 <r, where R 
and a are radius and surface density of the shell. Then, 
gravitational potential of this shell in a point with radius 
vector r 3 (from the center of the shell) is 



, . GnuRcl 
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r 3 



cosh 



2Tr(r 3 +R) 
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cosh 



27r(r 3 -i?) 



1 



(7) 



This formula demonstrates two features of the considered 
models. Firstly, we see that inside (r 3 < R) of the shell 
gravitational potential is not a constant. Thus, a test 
body undergoes an acceleration in contrast to the newto- 
nian case (and Birkhoff 's theorem of general relativity in 
four-dimensional spacetime which states that the metric 
inside an empty spherical cavity in the center of a spher- 
ically symmetric system is the Minkowski metric). Sec- 
ondly, this potential has a logarithmic divergency when 
r 3 - R: V (r 3 ) a - %2 [l - ^ In (2n\R - r 3 \/a)] 
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where we took into account R >> a and \R — r 3 \ << a. 
For example, in the case 2ttR — 10cm and 2tt\R — r 3 \ = 
10 -1 a, the deviation constitutes 2.3 ■ 10~ 4 and 2.3 ■ 10~ 18 
parts of the newtonian value —Giym/R for a — 10~ 3 cm 
and a = 10~ 17 cm, respectively. In principle, the for- 
mer estimate is not very small. However, it is very dif- 
ficult to set an experiment which satisfies the condition 
\R — r 3 \ << a. If the shell has a finite thickness, then 
the divergence disappears. 

Spherical shell The gravitational potential of a spher- 
ical shell of inner radius i?i and outer radius R 2 can be 
written as 



<p(r a ) 



Gn pa 
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cosh 
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(8) 

where p is a constant volume density of the shell: p = 
raj \^§-{R\ — Rf)]- It is useful to present this potential 
in the form of series. For example, inside (r 3 < Ri) of 
the shell it reads: 



<p{r 3 ) = 2G NP \ -tt (R 



( a \ I Zttk 




where we singled out zero mode k — which corresponds 
to the newtonian limit. It can be easily seen that this 
series does not diverge when r 3 — * R±. However, acceler- 
ation of a test body diverges when it approaches bound- 
aries Ri and R2 of the shell. It clearly follows e.g. from 
the form of the acceleration inside (r 3 < i?i) of the shell: 
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(10) 



and outside (r 3 > R2) of the shell: 

Gjym Gat pa 3 
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\ <<xj> ( —r 3 ) [h k (R 2 ) - h k (Rx)} < . 



(11) 



where h k {R) = ^Rcosh(^R) - sinh(^-R). Diver- 
gence of acceleration originates from the divergence of the 
series J2t^i an( ^ nas ^ ne f° rm ±2GAr/9aln Here, 
£ = \Ri,2 — r 3 \ and — ,+ corresponds to r 3 — > R\ and 
^3 — ► R2, respectively. In the case of a sphere [R\ = 
and i?2 = R) this divergence can be rewritten in the form 
■ ln(27r|i? — r 3 \/a). Similar to the case of the 
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" R? 2tvR 

infinitesimally thin shell, this deviation from the newto- 
nian acceleration — G^vm/i? 2 is also difficult to observe 



at experiments for considered above parameters. Eqs. 
(TIT)]) and (TTTj) show that acceleration changes the sign 
from negative outside of the shell to positive inside of 
the shell. This change happens within the shell. In the 
limit R\ — y 0, some of obtained above formulas can be 
applied for a solid sphere. 

Gravitational self-energy of a sphere For a sphere 
of constant volume density p, mass m and radius R, a 
gravitational self-energy is 



fe=i 



U = E/jvjl (12) 
-.r f a \ 3 V^ 1 \2nkR ( a \ 2 

where 



F k = 1 



27rfci?\ / 2nkR 

exp 

a J 



2irkR , (2nkR 
cosh 



sinh 



and Um = — 3G,/vra 2 /(5-R) is gravitational self-energy of 
a sphere in the newtonian limit @. In the case R» a, 
the difference AU = U — Un reads: 

„ 15 /a\ 2 2irt^ 1 rT 5 /a\ 2 

/c=i 

where we took into account X^feLi 1/& 2 = tt 2 /6. There- 
fore, this difference is suppressed by power law (with re- 
spect to ratio between a and R) but not exponentially as 
it is usually expected for Kaluza-Klein modes. Neverthe- 
less, for Sun with R w 7 ■ 10 10 cm this value is a negligible 
part of U N even if o a 10 _3 cm : AU w 4 • 10~ 29 C/Ar 
where U N & -2 ■ 10 48 erg @. 

Gravitational interaction of two spheres Energy of 
gravitational interaction between two spheres of constant 
volume densities p, p' , masses m, m and radiuses R,R' 
on a distance r 3 > R + R' reads 

u ^ = -^—^ 1 + 18 (^r) (2^) 



^^exp (-^^3) h k {R)h k (R')^ 



(13) 



The member of the series with k — 1 (first Kaluza-Klein 
mode) gives the main correction to the newtonian expres- 
sion and acquires the form of Yukawa potential. In this 
case, for the force of gravitational interaction between 
these two spheres we obtain: 



dU 
dr 3 



G^mm! 



9 / a \ 2 / a 
+ 2V2^R/ \2ttR 1 



2nr 3 



exp 



~-(r 3 -i?-i?') 
a 



(14) 
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where we made additional natural assumption R, R' > > 
a. If r — R — R' ss a, then we get an estimate: 

dU G N mm' ( 

- — f» g i 1 ( 15 ) 

+ — (sg)W 2 ?}- 

For example, in the case 27ri? = 2-7T-R' = 10cm and a sa 
10 _2 cm the correction is 1.68 • 10 -11 , which is difficult 
to observe at experiments. However, in the case of the 
internal space topology T d the correction term in (|T5|) 
acquires a prefactor s which satisfies the condition 1 < 
s < d and represents a number of extra dimensions with 
periods of the torus which are equal (or approximately 
equal) to a = maxa,. If all of a,i are equal to each other, 
s = d. Increasing the number of extra dimensions d, 
finally we arrive at the condition when the correction 
term becomes big enough to contradict with experimental 
data. Therefore, in this case we can get a limitation on 
a maximal number of extra dimensions for considered 
models. Certainly, the models with infinite number of 
extra dimensions with a = maxo, are forbidden. 

Conclusions We have considered generalization of the 
Newton's potential to the case of extra dimensions where 
multidimensional space has topology M D = R 3 x T d . It 
was shown that for model with one extra dimensions in 
the case of massive point source, the gravitational poten- 
tial can be expressed via compact and elegant formula 
([6]). This formula is valid for full range of variables r 3 



(magnitude of a radius vector in three dimensions) and £ 
(extra dimension) and has well known asymptotic behav- 
ior. Then, this formula was applied to an infinitesimally 
thin shell, a shell, a sphere and two spheres to get gravita- 
tional potentials and acceleration of a test body for these 
configurations and to compare obtained results with the 
known newtonian formulas. In some cases, obtained po- 
tentials and accelerations have logarithmic divergences 
near the boundaries of shells and spheres. Additionally, 
in contrast to the newtonian case, test bodies acceler- 
ate inside of shells. For each considered problems, we 
found deviations from the known newtonian expressions 
and show that for proposed parameters of the models it is 
difficult to observe these deviations at experiments. Nev- 
ertheless, if internal space has topology T d with approx- 
imately equal periods et^ of the torus in different dimen- 
sions, then we can get a limitation on maximal number of 
extra dimensions from experimental data with the help of 
formulas (fT^ . (fTS]) for gravitational interaction between 
two massive spheres. To conclude this letter, we want to 
stress that with the help of exact formulas of the form of 
(|6"|) we can solve multidimensional quantum Schrodinger 
equation. It opens a possibility to find a new gravitation- 
ally bound states: "dark matter atoms". We postpone 
this investigation to our forthcoming paper. 
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